Density of states of disordered Dirac particles: Infinitely many operators with 
negative scaling dimensions and freezing transitions 
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The global density of states (GDOS) close to the band center e = for a particle hopping on 
a square lattice and subjected to disorder that preserves the bipartite symmetry of the lattice 
is computed using field theoretical methods. The GDOS diverges like |e| _1 exp (— c| In e| K ) with 
k = 2/3 in agreement with a prediction by Motrunich et alt— and in disagreement with an older 
prediction by Gade- (ft = 1/2). 
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The problem of two-dimensional (2D) disordered Dirac 
particles, particles with a relativistic dispersion subjected 
to several types of quenched randomness, has attracted 
a lot of interest since 1994; 1 Indeed, it is related to a 
broad class of models such as the 2D random phase XY 
model, dirty d-wave superconductors, or the Chalker- 
Coddington model for the plateau transition in the in- 
teger quantum Hall effect, in which disorder plays an 
essential role. In spite of its relative simplicity, it can 
exhibit a rich variety of phenomena which appears to be 
specific to critical behavior induced by disorder. Multi- 
fractality is one example thereof. 

Multifractal scaling implies that the scaling regime in 
the vicinity of criticality is encoded by an infinite set of 
independent relevant scaling exponents, the multifractal 
spectrum in short. On general grounds, the multifractal 
spectrum must fulfill the condition of monotonicity. Im- 
plementing this condition is a challenge to the common 
lore of local field theoretical techniques. 

We consider a problem of Anderson localization de- 
fined on a 2D lattice where a particle is hopping on a 
square lattice in the background of 7r-flux phase and sub- 
jected to weak uncorrelated bond-disorder that preserves 
the sublattice symmetry and time-reversal invariance. 2 
The trademark of this model, the HWK model, is the 
sublattice symmetry present for any given realization of 
the disorder. Numerical studies offer strong evidences 
that typical wavefunctions at the band center are neither 
localized nor delocalized but are multifractal^ and that 
the density of states (DOS) is singular 3 . This disorder 
induced critical behavior at the band center is caused 
by the sublattice symmetry. After taking the continuum 
limit, the HWK model reduces to two flavors of Dirac 
particles subjected to an imaginary Abelian random vec- 
tor potential and a complex valued random mass. 

It is known that monotonicity of the multifractal spec- 



trum for the zero-modes of one flavor of 2D Dirac parti- 
cles subjected to a white-noise random vector potential 
is implemented in a rather dramatic fashion, namely by 
the termination of the multifractal spectrum^ In anal- 
ogy to thermodynamics, one can recast the termination 
of the multifractal spectrum as a freezing transition. 
From a field theoretical perspective, the termination of 
the multifractal spectrum in this model is closely re- 
lated to the existence of an infinite set of operators 
with negative scaling dimensions whose renormalization 
group (RG) flows are governed by a non-linear diffusion 
equation, the Kolmogorov-Petrovsky-Piscounov (KPP) 
equation&i In this article, we sketch how, within field 
theory, the global (i.e., self-averaging) density of states 
(GDOS) in the HWK model is governed by the same 
KPP equation and diverges at the band center e = 
as | er | 1 exp (— c| lne| K ) with k — 2/3 in agreement with a 
prediction by Motrunich et al. s and in disagreement with 



an older prediction by Gade 



1/2). 



We start from the continuum limit of the HWK model. 
The Hamiltonian is given by (ct's are the Pauli matrices 
and the unit 2x2 matrix) 
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i<r A + a 3 A 3 , 



(1) 



where A — [AijA^) £ R 2 is the purely imaginary ran- 
dom vector potential, Aq G M. the purely imaginary ran- 
dom scalar potential (chemical potential), and A 3 e R 
the random mass potential. We assume that all poten- 
tials are white-noise Gaussian distributed with vanishing 
means and with the variances specified by A„(x)A v (y) = 
g A S lll ,S(x - y), for /j,,u = 1,2, and A^(x)A v (y) = 
5(x - y), for /j,, v = 0, 3. Disorder averaging is de- 
noted by an overline. We employ the supersymmetric 
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(SUSY) method and introduce four-component Grass- 
mann fields ?/>, ^, and their bosonic partners j3, (3 to 
represent the resolvent (H HWK — e + iif)~~ of the HWK 
Hamiltonian in a given realization of the disorder. The 
advantage of the SUSY representation is that disorder 
averaging over the resolvent reduces to a Gaussian inte- 
gration and induces an effective interacting theory with 
action S £ -i n for the fields ip, ip, f3, and (3. 

Before characterizing the statistical distribution of the 
DOS, we need to understand the action S £ -iri=o- When 
the random mass perturbation is switched off, g M = 0, 
and at the band center, e — irj = 0, the resulting effec- 
tive interacting field theory S 1 * defines a conformal field 
theory (CFT) with gl(2\2) current algebra symmetry for 
any value of g A ~ It is also known that g M > is ex- 
actly marginal whereas g A is marginally relevant with 
the non-perturbative beta function (3 g ^ = (<? M £) 2 /(27r 2 ), 

C = (l + g M /7r) _1 pi£ More precisely, switching on a finite 
g M > demotes the CFT to the status of a nearly con- 
formal invariant theory (NCFT)Jfl Although conformal 
invariance is now broken, correlators for the currents gen- 
erating the gl (2 1 2) algebra can nevertheless be computed 
non-perturbatively in g A and g M . Furthermore, a weak 
form of Wick theorem still holds in the NCFT. This is 
all we will need to compute the GDOS. 

To calculate the energy dependence of the GDOS, 
which we assume to be identical to that of the disorder- 
average local density of states (LDOS) v m {e), one must 
reinstate a finite complex energy, i.e., a GL(2|2) sym- 
metry breaking perturbation An=i = ipip + P/3 is added 
to the action iS £ _i J?= o with the complex energy e — in 
as a coupling constant. The RG flow of the expectation 
value for tj){x)ij){x) + (3{x)(3{x) in the critical theory per- 
turbed by a finite but small energy under an infinitesimal 
rescaling a — * ae dl of the ultra-violet cutoff o gives the 
Callan-Symanzik equation obeyed by v av (e). 

The energy perturbation Ai is one member of an in- 
finite set of operators with negative anomalous scaling 
dimensions piA Other members, An — "mC^i) + "'> 
from this set are generated through the operator product 

N + N'\ 



x An> = 



expansion (OPE) 11 An > -~> \ — ( ^ 

. The anomalous scaling dimension x N of An is 



x N = ([l + 0(g M )]N 
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As the relevance of An grows with N or g A , all An 
must be taken into account on an equal footing when 
performing a RG analysis of i/ av (e) . 

The perturbed action needed to perform a consistent 
RG analysis of v m {e) is obtained by adding to S e -i V= o 
the perturbations -EL y iY (lI ''" 2 / (l2r ^( r )' An 
expansion of the one-point correlation function v a - v (e) up 
to second order in the "fugacities" Y N shows that it obeys 
a Callan-Symanzik equation characterized by infinitely 



many one-loop /3-functions j3 Y = dY N /dl, 
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(3 Yn = (2 - x N )Y N + tt ( m> ) Yn ' Y ' 
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(3) 



with the initial condition Y N (l = 0) = (e — ir))/a Xl ~ 2 , 
together with an anomalous scaling operator that cou- 
ples the expectation value of ip(x)'ip(x) + (3(x)(3(x) to all 
higher powers thereof. This Callan-Symanzik equation is 
thus an infinite set of coupled differential equations. 

To treat all An on equal footing, we define the gener- 
ating function 



Y n (l). 



(4) 



The RG equation for this generating function reduces to 
a non-linear diffusion equation of the KPP-type^ 
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dy)G + 2G(G - 1) 



(5) 



with the initial condition G(y, 0) = cxp (— ^ ^'^ e" 

A detailed study of the asymptotic [l — » oo) solution of 
Eq. (J5J suggests that the GDOS v aN {e) is governed by 
the "equivalent" Callan-Symanzik equation 11 
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The infinite set of coupling constants Y N together with 
the infinite dimensional anomalous scaling operator have 
been replaced by a single coupling constant z A , the dy- 
namical scaling exponent, that relates a change of the 
typical energy scale to a rescaling of length. From the 
perspective of Eq. |JSJ, z A is uniquely defined as the ve- 
locity of the asymptotic wave front solution, ft scales 
as g A for g A Q 2 < 2n whereas it scales as ■Jg~ A for 
g A (, 2 > 2tt. A freezing transition must necessarily take 
place at g A C 2 = 27r as g A is marginally relevant. Equa- 
tion © can be solved. The most interesting case is 
that of weak disorder g A (l — 0)£ 2 C 1 at the bare 
level. Sufficiently far away from the band center, i.e., 
when the sublattice symmetry is completely broken by 
the perturbing bare energy, the GDOS is flat. Upon 
approaching the band center, the sublattice symmetry 
manifests itself by a power-law behavior of the GDOS 
when the marginal relevance of the Abelian random vec- 
tor potential can be neglected. At lower energies the rel- 
evance of g A (l)( 2 ~ 1 < 2ir cannot be neglected anymore 
and the GDOS starts to increase in a Gade-like man- 
ner. Arbitrary close to the band center the relevance 
of g A (l)C 2 > 2-7T dominates everything else in that the 
random vector potential induces a strongly fluctuating 
LDOS (in a log-normal manner). In the strong disorder 
regime g A (l)C 2 > 27r, z a has undergone the freezing tran- 
sition that is responsible for zero-modes displaying local- 
ized behavior instead of multifractalityi and the GDOS 
diverging with the exponent k = 2/3 instead of k = 1/2&. 
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In summary, the HWK model is described in the con- 
tinuum limit by a NCFT in which we identified infinitely 
many operators with negative anomalous scaling dimen- 
sions. These operators reflect the broad distribution of 
the LDOS and enter in the field theoretical computa- 
tion of the GDOS through a Callan-Symanzik equation 
made up of infinitely many differential equations. From 
a detailed study of the RG equation obeyed by a gener- 



ating function for an infinite subset of relevant fugacitics 
we identified a scale dependent dynamical exponent that 
quantifies the change in the typical energy scale due to a 
rescaling of length. We then conjectured that the GDOS 
obeys an equivalent Callan-Symanzik equation made of a 
finite number of differential equation and whose solution 
agrees with the prediction from Motrunich et al. for the 
GDOS in the HWK model. 
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